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Generalisation of the Eyring-Kramers 
transition rate formula to irreversible 
diffusion processes 

Freddy Bouchet and Julien Reygner 


Abstract. In the small noise regime, the average transition time be¬ 
tween metastable states of a reversible diffusion process is described at 
the logarithmic scale by Arrhenius’ law. The Eyring-Kramers formula 
classically provides a subexponential prefactor to this large deviation 
estimate. For irreversible diffusion processes, the equivalent of Arrhe¬ 
nius’ law is given by the Freidlin-Wentzell theory. In this paper, we 
compute the associated prefactor and thereby generalise the Eyring- 
Kramers formula to irreversible diffusion processes. In our formula, the 
role of the potential is played by Freidlin-Wentzell’s quasipotential, and 
a correction depending on the non-Gibbsianness of the system along the 
instanton is highlighted. Our analysis relies on a WKB analysis of the 
quasistationary distribution of the process in metastable regions, and 
on a probabilistic study of the process in the neighbourhood of saddle- 
points of the quasipotential. 


1. Introduction 

1.1. The Eyring-Kramers formula 

Many equilibrium systems in statistical physics can be described by the over¬ 
damped diffusion of a particle in according to the solution (X^) t > 0 of the 
stochastic differential equation 

d X e t = -XU(Xl)dt + V2~e& W t (1.1) 

in M d , where ( Wt)t>o is a ci-dimensional Brownian motion, e > 0 is a tem¬ 
perature parameter and U : K. d —> R is the potential of the system. Upon 
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integrability assumptions on this potential, the process is reversible and er- 
godic with respect to the Gibbs measure 


/ 40 ) = J7 ex P 



Z £ := / exp 

iiCR 1 * 



( 1 . 2 ) 


which in particular implies that the particle visits infinitely often every region 
of the state space However, when e is small and U possesses several local 
minima, or potential wells, the particle typically remains stuck at the bottom 
of these wells over long times, and transitions between different wells are rare 
events of the system. This is the metastability phenomenon. 


AGi 


Figure 1 . Example of a double-well potential in dimension 
d = 2. The potential barrier to escape from the well corre¬ 
sponding to the local minimum X\ is Af7i = U(x+) — U{x i). 



It is of course of major interest to quantify the time scale over which 
such transitions between different metastable states occur. A first answer is 
provided by Arrhenius’ law [2j, asserting that the average exit time E[r £ ] 
from a potential well satisfies the logarithmic equivalence 


lime logE[r e ] = A U, (1.3) 

e4-0 


where AU is the minimal potential barrier that the particle has to climb in 
order to escape from the well, see Figure |T] In the sequel of the paper, we 
shall use the notation 

E|T ‘ l ,fo eXP (^) (L4) 

to refer to such a logarithmic equivalence. 

Let two potential wells, corresponding to local minima X\ and x^ , com¬ 
municate through a saddle-point (see Figure [Q. We make the generical 
assumption that the Hessian matrix Hess U (x \) is positive-definite, and that 
the Hessian matrix Hesst/(x*) has one negative eigenvalue — A+, and d — 1 
positive eigenvalues. Let us denote by the transition time between the 

two wells. Then the subexponential correction to Arrhenius’ law is given by 
the Eyring-Kramers formula 


| det Hess C/(x*)| 
det Hess U{x\) 


x exp 


(AU 




e 


(1.5) 
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The formula is named after Eyring and Kramers’ respective papers hsieii. 
Although we do not exactly know when this formula was first derived for 
dynamics in any dimensions, it already appeared for that case in Landauer 
and Swanson [32], and Langer [33] . 

The mathematical proof of this formula, as well as its extension to more 
general metastable situations, was first obtained by Bovier, Eckhoff, Gayrard, 
and Klein HU by means of a potential theoretic approach, and shortly after 
by Helffer, Klein and Nier m through Witten Laplacian analysis. We refer 
to the survey by Berglund [Bj for a review of mathematical approaches to the 
Eyring-Kramers formula, and highlight that the study of such formulas for 
general equilibrium systems remains an active field in mathematics, see for 
instance HUH El for degenerate potentials, [40] for the relation with func¬ 
tional inequalities, or [3 [EU H E2] for infinite-dimensional models described 
by stochastic partial differential equations and associated particle systems. 

1.2. Metastability out of equilibrium 

The purpose of this article is to establish a similar Eyring-Kramers formula 
for general diffusion processes 


d XI = b(Xf)dt + V2^a(XZ)dW t , 


( 1 . 6 ) 


that are not necessarily reversible, so that they can describe physical sys¬ 
tems out of equilibrium. In the small noise asymptotics, the particle typically 
fluctuates around the attractors of the zero-noise dynamics 


x t = b(x t ), 


(1.7) 


which are the metastable states of the system. In this paper, the zero-noise 
dynamics is called relaxation dynamics , and we shall assume that its only 
attractors are isolated points. We refer to [23] or [1(1, m n for a descrip¬ 
tion of the behaviour of more general nonequilibrium systems in the small 
noise regime. In our simple situation, the metastability of the process is ex¬ 
pressed by the fact that the particle spends most of its time fluctuating in the 
neighbourhood of attractors, where the law of its position is described by a 
quasistationary distribution, and rarely jumps from one attractor to another 
one, under the effect of the noise. The typical paths employed by the parti¬ 
cle for such fluctuations far from the typical behaviour are called fluctuation 
paths. 

At the logarithmic scale, the expected transition time between two at¬ 
tractors is described by the Freidlin-Wentzell theory [2D] , which extends Ar¬ 
rhenius’ law to irreversible processes. In this theory, the role of the potential 
in the reversible case is played by a quantity called quasipotential , constructed 
in terms of the action of fluctuation paths. In a simple situation where the 
relaxation dynamics (Oi possesses two attractors x\ and X 2 , the transition 
time from x\ to X 2 generally satisfies 



( 1 . 8 ) 
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where V{x\,-) is the quasipotential with respect to X\ and :c* is a relevant 
saddle-point of the quasipotential, located on the hypersurface separating the 
respective basins of attraction of x\ and X 2 ■ When e J, 0, the typical trajectory 
employed by the process to reach the saddle-point is deterministic. It is 
referred to as the instanton , or sometimes most probable escape path j38j , and 
is denoted by (pt)*€K- It connects X\ to cc* in infinite time, so that it satisfies 


(1.9) 


lim p t =x i, lim p t = x*. 



1.3. Main results and outline of the article 

The purpose of this article is to supply (11.81) with a subexponential prefactor, 
similarly to the Eyring-Kramers formula. We mostly rely on two ingredients: a 
perturbative calculation of the quasistationary distribution in the neighbour¬ 
hood of attractors, and a probabilistic study of the diffusion process (HI) in 
the neighbourhood of saddle-points. 

Let us explain our argument in further detail. In Section [2] we begin by 
recalling some notions from the Freicllin-Wentzell theory [JO] of large devia¬ 
tions for the diffusion process (HI). In particular, we give the definition of the 
quasipotential with respect to an attractor of the relaxation dynamics, and 
insist on the fact that this quasipotential can be identified when the vector 
field b admits a decomposition into a gradient part and a transverse force. 
We then discuss the formulation of this transverse decomposition in terms of 
a Hamilton-Jacobi equation for the quasipotential, which is related with the 
regularity of the quasipotential and plays a crucial role in the sequel of the 
paper. The study of the regularity of the quasipotential is a nontrivial math¬ 
ematical issue and leads to interesting physical phenomena [3]. Throughout 
the paper, we avoid considering such phenomena, and we shall rather estab¬ 
lish our main result under the assumption that the quasipotential be smooth 
in the neighbourhood of the instanton. Assuming that the quasipotential 
be smooth in a neighborhood of the instanton is not so restrictive, as this 
situation is thought of being generic, even in situations when Lagrangian 
singularities do exists in some part of the phase space. 

Sections [3] and 0] are dedicated to the study of the quasistationary dis¬ 
tribution of the diffusion process Ob in the neighbourhood of an attractor. 
In Section 01 we compute the prefactor to the ensemble measure, that is, the 
nonequilibrium stationary distribution of a particle diffusing in a vector field 
with a single global attractor. At the logarithmic scale, the Freidlin-Wentzell 
theory implies that the ensemble measure is equivalent to a Gibbs measure 
with the quasipotential playing the role of the potential. Using a classical 
WKB approach ,03 03 123 , 37 , [ 35 ], we show that the transport equation 
satisfied by the prefactor to this logarithmic equivalent can be solved us¬ 
ing the most probable fluctuation paths of the process as characteristics, or 
rays PH Hum], and thereby obtain an explicit expression for this prefactor 
which provides a correction taking into account the non-Gibbsianness of the 
system along these fluctuation paths. In Section 0 ] we fix a domain D cM. d 
attracted to a single equilibrium point of the relaxation dynamics H l and 
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perform a boundary layer analysis of the quasistationary distribution in D 
to compute the associated probability current at some point y £ dD. Inte¬ 
grating this current over dD yields the rate Aq St at which a particle escapes 
from D under the quasistationary distribution. 

In the bistable situation described in Subsection 11.21 the irreversible 
(or nonequilibrium) Eyring-Kramers formula is finally derived in Section [5] 
In contrast with Kramers’ method HU 133 [Ml iH HI], we do not apply a 
WKB approximation to the elliptic problem solved by the average transition 
time, but rather combine the results of Section 0] with a purely probabilistic 
study of the diffusion process in the neighbourhood of the saddle-point, in 
order to compute the rate at which a particle started from x\ reaches the 
neighbourhood of X 2 ■ Taking the inverse of this rate yields the following 
expression for the average transition time 




2-7T 


I det U* 


40 A* y det Hess x V{x\,x\) 
f V(xi,x 


exp 


✓»+oo 


F(p t )dt 


x exp 




( 1 . 10 ) 


In this irreversible Eyring-Kramers formula, H+ has to be understood as the 
suitable definition of the Hessian of the quasipotential at the saddle-point, 
so that the ratio of Hessian determinants is a straightforward generalisation 
of the equilibrium case, with the quasipotential playing the role of the po¬ 
tential. Besides, A+ is the positive eigenvalue of the Jacobian matrix D6(a;*), 
corresponding to the unstable direction of the vector field at the saddle-point. 
Note that, in general, A+ is not an eigenvalue of —if*, so that it depends on 
the whole vector field b and not only on its gradient part. The last correction 
with respect to the equilibrium Eyring-Kramers formula is the term involv¬ 
ing the integral of the function F along the instanton (pt)te R- This function 
F arises from the prefactor to the quasistationary distribution, and has to 
be interpreted as a measure of the non-Gibbsianness of the system, in the 
sense that F = 0 if and only if the Gibbs measure (L3 is stationary for the 
process ClU). 

Since our analysis is local, we expect this formula to extend to more 
general metastable situations, with more than two isolated attractors. 

Let us precise that similar nonequilibrium formulas have already been 
derived, although as far as we know, our probabilistic approach is new and 
leads to the original formula (11.101) . In specific restricted cases, Maier and 
Stein |37( 138) obtained a two-dimensional formula, while Ariel and Vanden- 
Eijnden PQ addressed the case of irreversible diffusion processes preserving the 
invariance of the Gibbs measure. Based on a purely WKB analysis, Schuss m 
obtained a preliminary version of (11.101) that does not explicitly involve the 
Hessian matrix if* of the quasipotential at the saddle-point. Besides, it seems 
to us that in none of these results, the nonequilibrium correction arising from 
the quasistationary behaviour was clearly identified in terms of the function 
F and of the instanton. 
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We finally underline that our analysis relies on the assumptions of local 
smoothness of the quasipotential around the instanton, and of nondegeneracy 
of the Hessian matrix of the quasipotential at the saddle-point. We believe 
that these assumptions are generic. For instance, for generic vector fields 
b, Lagrangian singularities of the quasipotential do not occur along the in¬ 
stanton. These assumptions essentially lead to the existence of a transverse 
decomposition for the vector field b , which somehow allows us to combine 
WKB approximations with notions from the theory of large deviations. In 
particular, the most probable fluctuation paths of the diffusion process are 
rather expressed in terms of this transverse decomposition than as the so¬ 
lution to the Hamilton equations in a 2d-dimensional phase space, which is 
the point of view of most previous papers D MM where singularities of the 
quasipotential are addressed. As we shall discuss below, our assumptions in 
particular exclude some specific phenomena, such as the grazing instantons 
observed in [181 SB] • 

1.4. Notations 

Throughout this paper, the scalar product in M. d is denoted by (■,■). Given 
a function / : ->• R, we write <%/ := ^§-/ and dijf := dx f dx . /• The 

Hessian matrix Hess/(a;) is defined by (Hess f(x))ij := dijf(x). Given a 
vector field h = (hi, ..., hd) : K. d K. d , the Jacobian matrix Dh(ir) is defined 
by (D h(x))ij := djfi{x). 


2. Quasipotential and transverse decomposition 

In this section, we introduce a few notions that are related with the Freidlin- 
Wentzell theory of the stochastic differential equation 

d XI = b(X e t )dt + VTea(Xl)dW t (2.1) 

in K. d , where b : R d —> and er : —> R dxm are smooth functions, and 

(Wt)t>o is a m-dimensional Brownian motion. This stochastic differential 
equation is written using the Ito convention. For all x € the diffusion 
matrix cr(a;)cr(a:) T G K. dxd is denoted by a(x). Throughout this section, it is 
assumed to be everywhere nondegenerate, in the sense that for all x G K. d , 
for all £ G \ {0}, (£, a(x)£) > 0. 

2.1. Generalities on the Freidlin-Wentzell theory 

We first introduce the Freidlin-Wentzell action functional, and the notion 
of quasipotential with respect to an equilibrium position of the relaxation 
dynamics. 

2.1.1. Relaxation dynamics and equilibrium positions. Given x G R rf , we 

will call relaxation dynamics started at x the path (ipf)t >o solution to the 
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(forward) Cauchy problem 

f tt = W)» * > o. , 2 2) 

\ r 0 = x- 

An equilibrium position of the relaxation dynamics is a point x £ W l 
such that b(x) = 0, so that iff = x for all t > 0. A stable equilibrium position 
is an equilibrium position x £ W l such that there exists an open subset 
O C R d containing x and such that, for all x £ O, iff converges to x when t 
grows to infinity. Note that this implies that the eigenvalues of the Jacobian 
matrix D b(x) have nonpositive real part. 

2.1.2. The Freidlin-Wentzell action functional. When the intensity e of the 

noise in ED vanishes, the sample-paths of the process )te[o,T] started at 
x £R d naturally converge to the relaxation trajectory (iff)te[o,T] > for all finite 
time T. The cornerstone of the Freidlin-Wentzell theory [SU] is the description 
of the fluctuations of the diffusion process around this deterministic limit 
through the action functional , generically defined for all Tj < T 2 and for all 
paths [r^Ta] by 

S[T u Tz],x{<f>) ■= T [ {4>t - b(<j>t), {(f>t - b{(f t )))dt (2.3) 

4 a t—Ti 

if cj) belongs to the set H 1 ([Ti,T 2 ]) of absolutely continuous paths such that 
(f> £ L 2 ([Ti,T 2 ]), and (j>r 1 = x; and 

S[T lt Ta],x(^) '■= +°° ( 2 - 4 ) 

otherwise. Then the Freidlin-Wentzell large deviation principle is the loga¬ 
rithmic equivalence 

VT > 0, P [Vf e [0, T],X e t ~ (j, t ] x exp . ( 2 . 5 ) 

We refer to [30) [TC| for a mathematical introduction to large deviation theory 
and for a rigorous formulation of this principle. 

2.1.3. Quasipotential. The convergence of the diffusion process to the relax¬ 
ation dynamics is naturally expressed by the fact that <S , [o,t],x(' ! /’ x ) vanishes 
for all T > 0. In fact, the action functional measures the difficulty for the 
diffusion process to deviate from its typical behaviour. This is best observed 
by introducing the quasipotential V (x, x) with respect to a stable equilibrium 
position x £ l d , defined for all x £ R d by 

V(x,x) := inf{Sp^.x (4>) ■ <hi =xAt 2 = a;,Ti < T 2 }. (2.6) 

Note that V(x,x) = 0, V(x,x) > 0 for all x £ R d , and the function x i->- 
V(x,x) is continuous on M d . 

2.2. Transverse decomposition of the drift 

In this subsection, we show that the quasipotential actually behaves like a 
potential for the irreversible dynamics. 
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2.2.1. Fluctuation dynamics and identification of the quasipotential. In this 
paragraph, we recall the framework of Jti] Section 4.3], where the quasipo¬ 
tential can be explicitly identified. Let us fix an open subset D of R d such 
that (b(y), n{y )) < 0 for all y £ dD , where n(y) refers to the outward normal 
vector to dD. The boundary is said to be noncharacteristic. This assumption 
ensures that the relaxation dynamics started at x £ D U dD cannot exit D. 

We furthermore assume that there exist smooth functions U : DUdD —> 
R and £ : D U dD —> R d such that: 

(i) for all x £ D, b{x) = —a(x)Vt/( x) + i{x) and (VU(x),£(x)} = 0, 

(ii) there exists x £ D such that, for all x £ D \ {5}, U(x) > U(x) and 
VC/ (x) ^ 0. 

The function U will be referred to as the potential. It can already be noted 
that, for all x £ D U dD, 

^ W) = -<a(V?)VE7(#), Vt/(#)} < 0, (2.7) 

so that U is a Lyapunov functional for the relaxation dynamics. As a conse¬ 
quence, when t grows to + 00 , ip x necessarily converges to x, which is thus a 
stable equilibrium position of the relaxation dynamics. 

Assuming this transverse decomposition of the vector field b, let us define 
the fluctuation dynamics terminated at x £ D U dD as the path (ipt)t< 0 
solving the (backward) Cauchy problem 

Rr = »Mwc%>n +t < 0, 

K=*. ( 1 

The fact that this path does not leave D is an assumption that has to be made 
on the shape of dD, for example, it is sufficient to assume that ( a(y)X7U(y ) + 
l{y),n{y)) > 0 for all y £ dD. When l = 0, the fluctuation dynamics is the 
time-reversal of the relaxation dynamics: 

Vt < 0, ^ t =V-v (2-9) 

In general, some properties related with time-reversal still hold true, for ex¬ 
ample the potential U remains a (backward) Lyapunov functional for the 
fluctuation dynamics, in the sense that 

> 0 . ( 2 . 10 ) 

This implies that, for all x £ D U dD, the fluctuation dynamics terminated 
at x satisfies 

lim ip* = x. (2-11) 

t—t—oo 

It then follows that (p x is an extremal path for (EH), in the sense that 
V(x,x) = U(x) — U{x) 

= \ f ^ b{v x t ),a{v x t )-\v x t - ))>dt, 

^ J t= — OO 


( 2 . 12 ) 
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see m Theorem 3.1, p. 100]. As a consequence, ip x describes the typical tra¬ 
jectory employed by the diffusion process to escape from the neighbourhood 
of x and reach x under the effect of a fluctuation of the noise — whence the 
designation fluctuation dynamics. 


2.2.2. Examples of transverse decompositions. Obviously, any gradient sys¬ 
tem, or overdamped Langevin process 

dX e t = — V£/(X t e )df + V2edW t , (2.13) 

provides a trivial situation where the drift vector field admits a transverse 
decomposition. In addition to modeling purposes in equilibrium statistical 
physics, such a stochastic differential equation has a computational interest 
as simulating its solution over long times allows one to sample from the 
associated Gibbs measure exp(— U/e) [33] . In this perspective, it is known 
that modifying the drift of the process in order to make it irreversible without 
changing its stationary measure generically improves its rate of convergence. 
A means to do so [251 US] El] is to consider the stochastic differential equation 

dX t e = -Vf7(X t e )df + JV[/(X t e )df + V^dW t , (2.14) 


with J being any constant antisymmetric matrix of size d. More generally, 
one can consider the stochastic differential equation 


d XI = -KXU{X^)dt + V2eadW t , 


(2.15) 


with the symmetric part \{K + K T ) of K given by a = <ra T . Then it is eas¬ 
ily checked that the Gibbs measure exp(— U/e) remains stationary. Besides, 
letting £(x) := \{K — K T )XU{x) yields the relation (XU(x), i(x)) = 0 every¬ 
where. This formulation includes the Kac-Zwanzig model and its generalisa¬ 
tion addressed by Ariel and Vanden-Eijnden [T] , as well as kinetic Langevin 
processes 


f d q t = p t dt, 

\ dp t = -X q u(q t )dt - ^p t dt + i/leydwt, 

for which U(q,p) — |p|/2 + u{q) and the diffusion matrix 

0 0 
0 7 


(2.16) 


(2.17) 


is degenerate. 

The examples above have the peculiarity of maintaining the Gibbs mea¬ 
sure exp (—U/e) stationary for the process, which is due to the fact that 
div i = 0 in these examples, but should not be the case for general irreversible 
processes. An instance of a vector field with transverse decomposition and 
for which the stationary distribution does not have an explicit expression 
is provided by the AB model from ff3], which displays similar features to 
infinite-dimensional models such as 2D Euler equations, the Vlasov equation, 
magneto-hydrodynamic equations, and the shallow-water equations. 
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2.2.3. Hamilton-Jacobi equation and genericity of the transverse decompo¬ 
sition. The relation 


Vz G D, (VU(x), t(x)) = 0 (2.18) 

of the transverse decomposition described above equivalently rewrites 

Mx G D, (VU{x),a(x)\7U(x)) + (b(x), VU(x)} = 0, (2.19) 

which is the Hamilton-Jacobi equation for the variational problem defining 
the quasipotential (1231) . As a consequence, even in the absence of an ex¬ 
plicit transverse decomposition, it can be proved that as soon as the func¬ 
tion x i—^ V(x,x) is C 1 , then it is a solution of this Hamilton-Jacobi equa¬ 
tion [20] Section 4.3]. Defining £(x) := b(x) + a(x)V x V(x, x), we deduce that 
{V x V(x, x),£(x)) = 0 for all x G R d . In other words, under some smoothness 
assumption on the quasipotential, the latter automatically provides an im¬ 
plicit transverse decomposition of the vector field b. Of course, this transverse 
decomposition is only useful in domains D satisfying the conditions of 42.2.11 
Singularities of the quasipotential, that is to say situations in which 
x K > V(x, x) fails to satisfy the Hamilton-Jacobi equation (12.191) . have been an 
important topic of research in nonequilibrium physics for several decades (241 
[251HBJ [371 [HI [ira, and have known a recent renewed interest due to their 
link with Lagrangian phase transitions in the macroscopic fluctuation theory 
of driven diffusive systems. We refer in particular to the recent review [J| 
by Baek and Kafri for a detailed survey. In this paper we shall not address 
this issue and generically assume that the quasipotential is C 1 in domains of 
interest. 


3. Prefactor for the stationary distribution 

This section is dedicated to the study of the stationary distribution of the 
diffusion process (A t e ) t >o solution to (12.11) . We shall therefore work under 
the standing assumption that for e > 0 small enough, this process possesses 
a unique stationary distribution P s £ t , and that this probability distribution 
possesses a smooth density p e st with respect to the Lebesgue measure on M. d . 
Practical conditions on the functions b and a ensuring this assumption can be 
found in [30] or [20]. Then the density p e st solves the stationary Fokker-Planck 
equation 

d d 

V.T G R d , 0 = e^ dij ( a ij (x)p € st (x )) - (, bi(x)p e st (x )). (3.1) 

i,j=1 i—1 


3.1. WKB approach 

When, for all x G a(x) is the identity matrix and b(x) = 

U(x) 


Z e := 


ceR d 


exp - 


—VU(x) with 
(3.2) 


da; < Too 
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it is well known that P s e t is the Gibbs measure defined by the density 


Plt( x ) = ex P 



(3.3) 


Assuming that U attains its global minimum at a unique point x G R d and 
that Hess U ( x ) is positive-definite, the Laplace approximation of Z e yields 
the equivalence 


Ptt(x) 


I det Hess U (x) 


ej.0 


( 27 re)° 


■ exp - 


U(x) — U{x) 


(3.4) 


In the general case, the WKB analysis of the stationary Fokker-Planck 
equation ED consists in looking for a solution of the form 


Pst 0*0 


C^x) 

e d/2 


exp 



(3.5) 


where the prefactor t (x)/e d ^ 2 is subexponential. This approach has been 
widely used |21 [391 ESUS] and it is known that injecting the ansatz above 
in (El and identifying the terms according to the powers of e yields the 
Hamilton-Jacobi equation (12.191) for U, and the transport equation 

(VC st , b + 2 aVU) + C st (div b + a : Hess U + 2 (A, Vt/}) = 0, (3.6) 


for the lowest order approximation of the prefactor 


C st (x) :=lim(7 s £ t (a:). (3.7) 

e^O 

In the transport equation (13.61) . the vector A{x) G R d is defined by 

d 

Ai{x) := J2 d j a ij{x)- (3-8) 

i =i 


3.2. Quasipotential in the eikonal equation 

In the context of WKB approximation, the Hamilton-Jacobi equation (12.191) 
for U is generally referred to as the eikonal equation. If the relaxation dy¬ 
namics possesses a unique stable equilibrium position x G such that 
x i ^ V(x,x) is C 1 on R d and V x V(x,x) ^ 0 for x ^ x, then the discus¬ 
sion of f )2.2.1l shows that V ( x , x) = U ( x ) — U(x) for all x G R. A probabilistic 
formulation of this result is the fact that the family of probability distri¬ 
butions {P s e t ;e > 0} satisfies a large deviation principle with rate function 
V(x,x) [21 Theorem 4.3 in Section 4.4], 


3.3. Explicit solution of the transport equation 

Under the assumptions of (12.2.11 let us denote l{x ) := b(x) + a(x)VU{x), 
with Vl/(i) = S7 x V(x,x). Then the transport equation (13.61) rewrites 

(VUst, aWU + £) + C st (div t + (A, Vt/)) = 0. (3.9) 

In the WKB approach, this equation is usually solved by the method 
of characteristics, or rays ED ESI- The main interest of the formulation of 
the problem in terms of large deviation theory is that these rays are exactly 
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given by the trajectories of the fluctuation dynamics. Indeed, the equation 
above rewrites 

(V log C s t, aVU + £) = —F, (3.10) 

where 

F( x) := div £{x) + (A(x),VU(x)}. (3.11) 

We note that if F = 0, then the Gibbs measure defined by (EH1) is stationary 
for the process, although the latter need not be reversible. In other words, F 
is a measure of non-Gibbsianness of the system. 

Recalling that a VIZ + £ is the vector field driving the fluctuation dy¬ 
namics, we obtain that, for all x £ R d , for all t < 0, 

A log C s t(tpt) = (V log C s t(<Pt), + %?)> = -F{tf), (3.12) 

so that integrating this identity for t £ (—oo, 0] finally yields 

C st (x) = C st (x) exp ^ F(<pt)dtj . (3.13) 


Under the assumption that HessU(x) = Hess x V(x, x) be positive-definite, 
performing a Laplace approximation in the normalisation condition for p e st 
yields 


Cst (5) 


I det Hess U (x) 


(27r)° 


(3.14) 


and we finally come up with the formula 


Pstfa) 


I det Hess x U(a:,a;) 


e4.0 


(27re) Q 


exp - 


V(x, x) 


F(tf)dt (3.15) 


for the stationary density. It is striking that it only differs from the for¬ 
mula ( 13 . 41 ) for the potential case through the integral term along the fluctu¬ 
ation dynamics, which therefore has to be interpreted as the accumulation 
along the fluctuation dynamics of some non-Gibbsianness of the system, mea¬ 
sured by the function F. 


4. Quasistationary exit rate of a domain 

In this section, we consider an open subset D C R d with noncharacteristic 
boundary dD and such that all the trajectories of the relaxation dynamics 
started in DUdD converge to a unique stable equilibrium position x £ D. As 
a consequence, when e is small, the diffusion process X e typically fluctuates 
around the point x, and this process hitting dD is a rare event. To quantify 
its time scale, let us introduce 

T dD '■= inf {£ > 0 ,X t e £ dD}. 


(4.1) 
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The Freidlin-Wentzell theory [2Dj Section 4.4] asserts that Tq D is of order 
exp (V(x,dD)/e), where V(x,dD) := min ye QuV{x,y). We therefore define 
the quasistationary phase as the range of times 


1 < t <C exp 




(4.2) 


during which the process forgets its initial position but typically remains 
stuck in a neighbourhood of x. 

The purpose of this section is to obtain a formula for the exit rate from D 
during the quasistationary phase. This exit rate is defined in Subsection l4.ll 
It relies on the solution to the stationary Fokker-Planck equation in D with 
absorbing boundary condition, which is well approximated in the range of 
times (14.21) by the quasistationary distribution introduced in Subsection 14.21 
The exit rate is computed through a boundary layer approximation for the 
quasistationary distribution in Subsection 14.31 


4.1. Probability current and exit rate 

Assume that the particle is absorbed, or killed, when it reaches the boundary 
of D. Then the law of its position, defined by 

Vf >0, Vse D, pl hs (t, x) \= P [X e t ~ x , tq D > t] , (4.3) 

satisfies the Fokker-Planck equation 
d d 

dtPlbs = e di ’3 ( a bPabs) - J2 di ( 6 *Pabs). t > xeD, (4.4) 
i,j= 1 i= 1 

supplemented with the absorbing boundary condition 

PabsO^) = 0, xedD, (4.5) 

see [22 Section 5.3.2], 

There is of course a loss of probability as time evolves, in the sense that 

/ Pabs(*> = P I T dD > *] ( 4 -6) 

J x£D 

decreases, corresponding to the probability that the particle exits D. This loss 
is measured by the probability current j £ (t, x), defined by the formulation of 
the Fokker-Planck equation (|4.4I) as the conservation law 

dtPl hs (t,x) + div j e (t,x ) = 0, (4.7) 

so that 

fit, x) = -e (A{x)pl hs + a{x)Xpi hs ) + b{x)pi hs , (4.8) 

for allt > 0 and x £ D. Then, the rate of exit from D is given by the flux of 
j e through dD , namely 

A e (i) := [ ( f(t,y),n(y))dS{y ), (4.9) 

J yEdD 

where d S(y) is the surface element and n{y) is the outward normal to dD , see 
for instance [3SJ Sections 2.5 and 2.6], [^] Section 10.2] or (2T] Section 5.4], 
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4.2. Quasistationary distribution 

Let us define the quasistationary distribution p qst associated with D by the 
so-called Yaglom limit of the conditional probability 



so that for t > 1, p\ hs {t, x) ~ p e qst (x)P [t% d > t}. On the other hand, on the 
time scale (14.21) . P [r| D > t] remains close to 1. As a consequence, during 
the quasistationary phase, the distribution p' bs (t, x) is well approximated by 
the quasistationary distribution p qst (x), and in particular it does not vary in 
time. 

During this phase, the exit rate thus writes 



(4.11) 


where the quasistationary current j qst is defined by approximating p| bs (f, x) 
with p qst (x) in (14.81) . On the boundary of D , the absorbing condition reduces 
its expression to 


Vy&dD, Jqst(y) = —ea(l/)Vpq St (y). 


(4.12) 


We compute this quasistationary current in the next subsection. 

4.3. Boundary layer approximation and computation of the quasistationary 
current 

The purpose of this paragraph is to compute the quasistationary current j qst 
on dD , in order to derive a formula for Aq St . 

4.3.1. Assumptions on the domain. We assume that the quasipotential x h-> 
V(x,x) is C 1 , and that, in addition to the fact that the boundary dD be 
noncharacteristic, it satisfies the condition that for all y € dD , 


i~a(y)V x V(x,y) + e(y),n(y)) < 0 and {\7 x V(x,y),n{y)) > 0, (4.13) 


where we have defined £{x) := b(x) + a(x)'V X V(x , x) and n{y) is the outward 
normal vector at y £ dD. This ensures that the results of T2.2.1I hold with 
U(x) = V(x,x), up to an additive constant. 

4.3.2. Boundary layer approximation. Conditioning the sample-paths of X e 
with respect to the event \tq D > t} in the definition of the quasistation¬ 
ary distribution induces a repulsive effect of the boundary on the particle. 
When the latter is far from dD, it does not feel this effect and therefore the 
quasistationary distribution in the bulk of D can be approximated by the 
expression (13.151) of the stationary distribution of a particle diffusing in a 
single-well potential given by the restriction of U to D. In the sequel, this 
distribution is referred to as the ensemble measure 



(4.14) 


where we recall that C s t{x) is defined by (13.131) and (13.141) . 
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On the other hand, the repulsion effect makes the quasistationary den¬ 
sity vanish in the neighbourhood of dD , where we therefore have to perform 
a boundary layer approximation to obtain the expression of p^ st {32] EHl IHJ ■ 
The thickness of this layer is determined by the fact that boundary layer ef¬ 
fects must be caused by a subexponential prefactor; in other words, if x £ D 
is close to dD and y e dD denotes its orthogonal projection on dD , then x 
is out of the boundary layer as soon as exp(— V(x,x)/e) exp(—V(x,y)/e). 

Writing, for \x — y\ <C 1, 

V (, x, x) ~ V(x, y) + (’ V X V(x, y),x-y) 

~ V(x,y) - \x — y\(V x V(x,y),n(y)), 

and using the assumption (14.131) that (' V x V(x,y),n(y )} > 0, we deduce that 
exp(— V(x,x)/e) exp(— V(x, y)/e) as soon as \x — y\ e, so that the 
boundary layer develops on a characteristic length scale e. 

This discussion leads us to assume that in the vicinity of dD, the qua¬ 
sistationary distribution writes 

riuw=(-bfe>), (4.16) 

where n d(x) is the orthogonal projection of x onto dD, rj{x) is the distance 
between x and ttd(x). Matching the quasistationary distribution with the 
ensemble measure p| ns in the bulk of D yields the boundary condition 

Vy e dD, lim C h \{y, r) = C st {y). (4-17) 

r —»+oo 

We point out the fact that the ansatz (14.161) is well suited to determine 
the leading order approximation of the quasistationary distribution, which is 
enough for our purpose. If one wants to obtain corrections at next orders, 
one should rather assume the prefactor to be a function of both x and i](x). 
This would lead to compatibility conditions to be solved. 

4.3.3. Computation of jq St (y). Differentiating the expression (14.161) with re¬ 
spect to x and injecting the result into the expression (14.121) of the quasista¬ 
tionary current on dD, we obtain 

jqstiv) - ex P ( 4 - 18 ) 

at the lowest order in e. Identifying with —n{y), we conclude that 

jqstiv) = ^a(yMy)^r(y,o)exp ^ ) ■ ( 4 - 19 ) 

We now compute ^§}r(y, 0). In this purpose, we fix y £ dD, r > 0 and 
let x := y — ern(y) € D. Then injecting the expression (14.161) of Pq St (x) into 
the Fokker-Planck equation (031) yields, at the lowest order in e, 

(n(y),a(y)n(y))^^-(y,r) + (a(y)V x V(x,y) + £(y),n(y))^^-(y,r) = 0. 

(4.20) 
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Note that the nondegeneracy assumption on a implies that (n(y ), a{y)n{y)) > 
0. Let 


viv) ■= 


(a(y)V x V(x,y) + £(y),n(y)) 


(n(y),a(y)n(y)) 

By (Id. Kill . fi(y) > 0. As a consequence, the equation (Id. 201) integrates as 


(4.21) 


f)n u , 1 — p-v{y)r 

Cb\(y,r) = Cbi(y,0) + —— (y,0) - — -, (4.22) 

or y{y) 

and the boundary conditions Cu(y, 0) = 0, which follows from (14.51) . as well 
as (14.171) imply 

^r(j/>0) = y{y)C st (y). (4.23) 

We deduce that the quasistationary current writes 


j qst (y ) 


Cst(y) 
f d /2 


exp 


V(x,x)\ ( a(y)\y x V(x, y) + l(y), n{y)) 
{n(y),a(y)n(y)) 


a(y)n(y). 

(4.24) 


4.3.4. Conclusion. Injecting the expression of the quasistationary current 
into the definition (14.111) of the quasistationary exit rate yields 

Aq St = f Cst }/ J 2 ex P (~ {a{y)^xV{x,y) + £(y),n(y))dS(y). 

Jy£dD ^ \ ^ J 

(4.25) 

One recognises here the expression of the ensemble measure Pg ns on the one 
hand, and of the vector field aV X V + t driving the fluctuation dynamics on 
the other hand, so that the quasistationary exit rate rewrites 

A qst = / (^t\t=o,n{y))pl ns (y)dS(y). (4.26) 

JyGdD 

This brings forth the final interpretation that the escape of the particle from 
D is governed by the flux through dD of the fluctuation dynamics, weighted 
by the ensemble measure Pg ns of the particle diffusing in D. 


5. The irreversible Eyring-Kramers formula 

We finally address the situation where the process X e is metastable in the 
neighbourhood of isolated attractors, and compute the average transition 
time between these attractors. As is discussed in the introduction of this ar¬ 
ticle, our analysis is local and can therefore be reduced to the simple situation 
of a bistable system. We therefore assume that the relaxation dynamics pos¬ 
sesses two equilibrium points X\ and x%, the respective basins of attraction of 
which are separated by a smooth hypersurface S containing a unique saddle- 
point a;*. As is underlined in [3B], the separating surface S is characteristic in 
the sense that (b(y), n(y)) = 0 for all y G S, which makes the boundary layer 
approach developed in the previous section more delicate, see also [um]. W e 
will therefore not follow this approach directly, and rather proceed as follows. 
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Consider a particle with initial position x\ and evolving according to 
the diffusion process m • We first use the Freidlin-Wentzell theory to show 
that, when e is small, the point at which the particle hits S is concentrated 
around 2*. This allows us to look at the linearised version of the dynamics 
in the neighbourhood of the saddle-point. For this dynamics, the separating 
surface between the two attractors is the tangent hyperplane to S at 2*, 
which we denote by So- For rj <C 1 , we denote by S v the hyperplane parallel 
to S 0 and located at a distance of order 77 of the latter, in the direction of 21, 
so that a particle starting from 21 needs to cross S v before reaching Sq. We 
then compute independently: 

(i) the rate at which particles arrive on S r}1 

(ii) the probability that a particle started from S^ drifts away from the 
saddle-point toward 22 instead of 21. 

Matching these two quantities yields the rate at which a particle started 
from 21 goes to the neighbourhood of 22, which is the inverse of the average 
transition time we are willing to compute. 

5.1. Bistability and instanton 

We first describe the bistability of the system and introduce a few notations. 

5.1.1. Vector field b. The relaxation dynamics is assumed to possess exactly 
two stable equilibrium positions aq and 22, with respective basins of attrac¬ 
tion D\ and D 2 separated by a hypersurface S = dD 1 = dD 2 . This hypersur¬ 
face is stable for the relaxation dynamics, and we assume that, for all 2 G S, 
ipt converges to the equilibrium position 2* G S , which is therefore stable in 
d — 1 directions and unstable in the remaining direction. See Figure [21 


S 


Figure 2 . Flow lines of the relaxation dynamics for a 
bistable system. 

We first prove that 

V(2i,2*) = inf V(xi,y), ( 5 . 1 ) 

y&S 

so that 2* behaves like a saddle-point for the quasipotential. To this aim, we 
fix y £ S and construct a continuous path (<(> t ) tg R as follows: 

• the part (ct>t)t <0 is an arbitrary continuous path joining aq to y, 

• the part (<l>t)t>o is given by the relaxation dynamics of from y to 2*. 




18 


Freddy Bouchet and Julien Reygner 


Then the part corresponding to the relaxation dynamics does not contribute 
to the action of (cf>t)te R> so that taking the infimum over all such paths we 
deduce that V{x\,y) > V(xi,x*). 

As a consequence, when e is small, the diffusion process typically fluc¬ 
tuates in the neighbourhood of a local equilibrium for long times, and then 
jumps to the neighbourhood of the other local equilibrium. If the quasipo¬ 
tential reaches its minimum on S only at the saddle-point, then the typical 
trajectory employed for this transition crosses S in the neighbourhood of x*. 


5.1.2. Geometry of the dynamics around the saddle-point. Let us denote 
M * := D5(x*) the Jacobian matrix of b at x*. On account of the assumptions 
on the vector field b around x*, the matrix M* has d — 1 eigenvalues with 
negative real part, and one nonnegative real eigenvalue. The corresponding 
stable spaces satisfy the following properties. 

(i) There exists a hyperplane V- such that M* V- C V- and, for all x € V-, 
exp(tM*)x converges to 0 when t —> +oo. The tangent hyperplane So 
to S at x* writes So = x* + VL. 

(ii) The nonnegative real eigenvalue is denoted by Aij_, it is assumed to 
be positive. The corresponding unit eigenvector is denoted by v+, we 
prescribe it to be oriented toward Lfy 

The normal vector to VL pointing in the direction of Zfy is denoted n*, its 
angle with u + is denoted 6. Given x £ M d , let us define £+(x) as the unique 
real number such that x — x* — £ + (x)v+ £ VL. We note that, for all x £ M. d , 


<+ 0*0 


(x - x*, n*) 
cosf? 


These notations are summarised on Figure [3] 


(5.2) 


5.1.3. Saddle-point and instanton. The instanton is the typical trajectory 
employed by the diffusion process to reach the saddle-point starting from xi. 
In general, it takes an infinite amount of time to leave xi, and an infinite 
amount of time to reach x*, so that the instanton has to be though of as a 
trajectory (p t )te R such that 

lim p t = x i, lim p t = x*, (5.3) 

t—>— oo £—)•+ oo 

for which the action 

1 f + °° 

j / (Pt ~ HPt),a(Pt)~ (Pt ~ b(p t )))dt (5.4) 

^ J t=— OO 

is minimal and worth R(xi,x*). In the sequel, we shall make the following 
assumptions: 

(i) there is a unique such path (up to time translations); 

(ii) for all t £ K, the quasipotential V (xi, •) is smooth in the neighbourhood 
of pt, so that it satisfies the Hamilton-Jacobi relation 

(V x V{x u x),V x V{x u x)) + (b{x),V x V{x i,x)) = 0; 


(5.5) 
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(iii) when t —> +oo, Hess x V{xi 1 pt) converges to some symmetric nonde¬ 
generate matrix H+ with d — 1 positive eigenvalues and one negative 
eigenvalue. 

Defining U{x) := V(xi,x) and £(x) := b(x) + a{x)S7U{x) in the neighbour¬ 
hood of each p tl we deduce from the same arguments as in Section [5] that 
{pt)t>o is a heterocline trajectory of the dynamical system 

p t = a(pt)S7U(p t ) + (5.6) 

connecting the equilibrium positions X\ and x*. In particular, for all t £ R, 
the fluctuation trajectory terminated at p t coincides with the instanton in 
the sense that 

Vs < 0, p* = p s+t . (5.7) 

Differentiating the transverse relation (VU(x), i(x)) =0 twice and eval¬ 
uating the result in p t for t —> +oo yields 

+ DjH+ = 0, (5.8) 

where D* := lim^+oo D £(pt) satisfies 

Ai* = + D*, a* := a(x*). (5.9) 

We finally note for further purpose that multiplying (15.81) by ff * _1 on both 
sides yields 

+ H ir ~ 1 D ir r = 0. (5.10) 


5.1.4. Linearisation of the fluctuation dynamics. The linearisation of the fluc¬ 
tuation dynamics around the saddle-point writes 


x = iV*(x — x*), TV* := + D*. (5.11) 


The flow lines of this dynamics are plotted in blue on Figure El In order to 
provide a justification of this picture, let us detail explicit computations in 
dimension d = 2, with a* being the identity matrix. 

We first change the coordinates so that x* = 0 and 




Pi 0 \ 
0 P2 ) 


(5.12) 


with pi < 0 < p 2 - Then (15751) implies that there exists afE such that 




0 a \ 

ap 0 J 


(5.13) 


with p \= —pi/p 2 > 0. Let us denote p := p 2 > 0, so that pi = —pp. A 
direct computation shows that the eigenvalues of M* are 


A* + := 


~p( i — p) + m(i + P)\J 1 + 


Apex 2 

M 2 (l +p ) 2 


AC := 


“M 1 ~ P) ~ M 1 + P)\l l + + p yi 


> PP > o, 


< — p < 0 , 


(5.14) 


and we denote by v+ and V- the respective associated eigenvectors. 
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Clearly, the eigenvalues of TV* are — and — A* , and we denote by v' + 
and v'_ the respective associated eigenvectors. The direction of the instanton 
incoming at the saddle-point is given by the eigenvector v' + associated with 
the negative eigenvalue —A* of TV*. A tedious but straightforward computa¬ 
tion shows that the angle 7 between v' + and V- satisfies 


sin( 7 ) = 


1 

y/l + (a/^ 


(5.15) 


which leads to the first conclusion that the assumption that H * be nonde¬ 
generate prevents the incoming direction of the instanton at the saddle-point 
to be parallel to Sq. In contrast with Maier and Stein’s results [35], this con¬ 
dition does not depend on the value of the ratio A+/|A* |, but rather on the 
sole assumption that \x > 0 . 

Likewise, the angle 7 ' between V- and v'_ vanishes if and only if a = 0, 
which leads to the second conclusion that if 7A* yT 0, then a region of D\ does 
not contain any fluctuation trajectory emanating from Xi: this is the clas¬ 
sically forbidden wedge evidenced by Maier and Stein [55], which is hatched 
on Figure [3] In this area, the quasipotential V{x\,-) is not described by a 
transverse decomposition. Therefore the function y > V[x \, y ) is generically 
not smooth on the hyperplane So, which is the reason why the distribution of 
the point at which the Brownian particle reaches S first is not Gaussian but 
rather skewed [55]. More generally, it makes Kramers’ analysis of the proba¬ 
bility current across S more delicate than the study presented in Section 0] 
which is the reason why we shall rather work on the hyperplane S v defined 
in (15.1611 below. 


5.2. Sketch of the argument 

We are interested in the probability that, during the quasistationary phase 
in D 1 , the particle crosses the separating hyperplane So and drifts toward 
X 2 . Following the approach described in the introduction of this section, we 
decompose this event into two parts, and compute separately: 

(i) the probability that the particle first reaches the affine hyperplane 

S v ■= x* - VV+ + V-= {C+(z) = — 77 }, (5.16) 

where 77 -C 1 will be specified below (see Figure 0]), 

(ii) the probability that a particle started from some point y £ S v actually 
crosses Sq and leaves toward X 2 - 

We want to apply the results of Section 0] in order to address the first 
point, and take as a domain D v the set {C+(a ; ) < — 77 }. Although the as¬ 
sumption (14.131) of Section 0] is not satisfied everywhere on S v , it still holds 
in the neighbourhood of the point y £ S v at which the instanton intersects 
S v . This is sufficient for our analysis, since the probability current essentially 
concentrates around this point. Then using the assumption of smoothness 
of V{x,y) in the neighbourhood of y, we replace V(xi,y) with its harmonic 
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Figure 3. The stable manifold Sq = x* + V- and the un¬ 
stable direction u+ of the linearised relaxation dynamics are 
plotted with thick red lines. Thin red lines represent flow 
lines of the linearised dynamics x = M* (x — x*). On this 
example, the coordinate £+(x) is negative. The thick blue 
line is the incoming direction v' + of the instanton. Thin blue 
lines represent flow lines of the linearised fluctuation dynam¬ 
ics. The hatched area is the classically forbidden wedge, in 
which the linearised fluctuation dynamics does not describe 
fluctuation trajectories emanating from xj.. 


approximation 

V(xi,y) ~ E(xi,x*) + - x*,#*(y - x*)), (5.17) 

where we underline that the distance \y— x*| is of order 77 -C 1. We then deduce 
from (14.241) that the probability that the particle crosses the hyperplane at 
some point y £ S v is given by 

,. e , , . Cst(x*) ( V(xi,x*) (y-x*,if*( 2 /-x*)) 

Oqst (V),n*) = ed/2 ex P --- Ye - 

x (IV*(y - x*),n*). 

As a consequence, for v / e <C 77, this distribution is approximately Gaussian, 
see Figure [I] The prefactor is computed in Subsection 15.41 below. 
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Figure 4. The background gray lines are the equipotential 
lines in the neighbourhood of the saddle-point for the qua¬ 
dratic potential ^(y — x+,H+(y — a;*)). The hyperplane S v 
introduced in Subsection 15.21 is plotted with the thick blue 
line. It is tangent to the equipotential line at the point y, 
determined in T5.4.1I The thin blue line is the instanton tra¬ 
jectory, it intersects S v at y. The distribution of the current 
on Sr) is plotted in light blue, it is approximately Gaussian 
around y and has a standard deviation of the order of sft. 


The second point is addressed by investigating the behaviour of the 
linearised diffusion process 

dXl = M*(X t e - x*)dt + V2ea{x±)dW t (5.19) 

far from the saddle-point. The particle following this linearised process started 
at a point of distance of order rj of is more likely to drift away in the di¬ 
rection of xi if its coordinate <(_|_(X t e ) hits —Z before hitting Z for 2 > 77 , 
see Figure [3] We therefore define the commitor function of the linearised 
dynamics by 

<?%):= lim P[f|<fl z ], (5.20) 

z^y+oo 

where Xq = y , and 

■■= inf{* > 0 : C+(*t) = ±Z}. 


(5.21) 
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The commitor function is computed in Subsection 15.31 



Figure 5. Typical and non-typical behaviour of the lin¬ 
earised diffusion process: when the particle starts on S v , the 
drifts tends to bring it toward negative coordinates (as mea¬ 
sured by £+), so that most trajectories hit {£+ = — Z} be¬ 
fore hitting {£_)_ = Z}. The commitor function describes the 
probability for a particle to perform the rare event of hitting 
{C+ = Z} before {C+ = -Z}. 


Combining the definition of q e (y) with (15.181) . we may finally conclude 
that the probability that the particle crosses the saddle-point during the 
quasistationary phase is given by 


C s t(a;*) 

fd /2 


■ exp - 


V(xi,Xi,) 


/ <f{y ) exp ( - 

J y £ Sri V 


{■y - x±,H+{y - X*)) 
2e 


(5.22) 


(N*(y ~ ir*),n*)d S(y), 


for rj satisfying y/e -C y <C 1. Then the average transition time is the inverse 
of this rate. 


5.3. Computation of the commitor function 

Let us fix y € and Z > 0 with Z > |£ + (y)|. For all t> 0, let us define 

4 ■■= (5.23) 
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where (A'|)t>o refers to the linearised dynamics started at y. Let us show 
that {zf ) t >o solves a one-dimensional stochastic differential equation, so that 
P[t z < ft z \ can be explicitly computed. On the one hand, using the decom¬ 
position of Xf — x * on v+ and V -, we obtain 


C+(M*(Xf — a:*)) = (5.24) 

while on the other hand, shows that the process (^ + (a*W t )) t >o is a real¬ 
valued Brownian motion with variance (a*n*, n*}/cos 2 0. As a consequence, 
{zl)t >o is the (repulsive) Ornstein-Uhlenbeck process 


dw t 


dzt = \* + ztdt + 

* + £ cos e 

and for all z £ [— Z, Z], the quantity 

q e z (z) :=P[t| < fl z \z e 0 = z] 
is obtained by solving the elliptic problem 

= »■ ^ I- 2 - 2 '’ 

with the boundary conditions 

= 0 , 

Integrating this equation yields 


Qz( z ) = 


L 


exp 


Z=-z 


f z (Z) = 1. 




2 e(a*n*,n*) 


A1 cos 2 9 „ 

exp 1 


2 e(a*n*,n*) 


dC 


(5.25) 

(5.26) 

(5.27) 

(5.28) 

(5.29) 


J £——Z 

and taking the limit when Z grows to infinity leads to the expression 

rC+(v) 


q e (y) = 


Afj cos 2 9 


27re(a*n*,n*) J i= _ c 
of the commitor function. 


exp - 


A+ cos 2 9 
2 e(a*n*,n*)' 




(5.30) 


5.4. Integration over the hyperplane 

The results of the previous subsection show that the commitor function q e (y) 
is uniform over y £ so that it can be taken out of the integral in the 
right-hand side of (15.221) . It therefore remains to compute 

I e v -= J s ex P (AT*(j/-a;*),n*)dS(j/). (5.31) 

We proceed in the following three steps. In 45.4.11 we determine the point 
at which the exponential term assumes its maximum, and highlight the link 
with the instanton trajectory. In 45.4.2l we express the value of I* in terms of 
the determinant of the restriction of the quadratic form v <—> (v. H*v) to V-. 
In 45.4.31 we connect this determinant with the determinant of H * and the 
eigenvalue A+. 
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5.4.1. Minimum of the potential along S v . We let y be the solution to the 
minimisation problem 

mi p -£*))> (5.32) 

yeSr, Z 

which rewrites 

min -(-i]v + +V- ) H ie (-i]v + + V-)), (5.33) 

(v- ,n*)=0 Z 

after the change of variable y = x* — r]v+ + V— . Defining the Lagrangian 

C{v~, A) := ^{-T]v+ + V-,H+(-r]v + +V-)) + A(u_,n*}, (5.34) 

the optimality condition 

V v _C = 0 (5.35) 

implies that there exists A £ K. such that H+(y — x*) + An* = 0, that is to 
say 

y = x* - A17* _1 n*. (5.36) 

Let us prove that A = —A +77 cos 0/(a*n*, n*}, and that y — x* is an 
eigenvector of the matrix IV* for the eigenvalue —At. We first proceed to 
prove that 

IV*Lf* _ 1 n* = —A ;F*- 1 n*. (5.37) 

To this aim, we write 

IV*iJ* _ 1 n* = (a*tZ* + -D*)iJ* _ 1 n* = a*n* - iJ* _ 1 D* T n* = -U* _ 1 M* T n*, 

(5.38) 

where we have used (15.101) for the second equality. Now for all u £ M d , using 
the decomposition of u on u_|_ and VL yields 

(M* T n*, u) = (n*, A+C + (it)i> + ) = A+(n*, u), (5.39) 

thanks to (IQl) . As a consequence, 

M* T n* = A * n*, (5.40) 

which completes the proof of (15.371) . We now turn to the evaluation of A, and 
first deduce from the series of equalities 

(y - x*, i/*(y - x*)) = -A(y - x*, n*) (5.41) 

and 

(j/-x*,fl*(j/-x*)) = --— (AT*(y-x*),fl'*(y-x*)) 

A + 

= - —(a*77*(y-x*),i7*(y-x*)) ( 5 . 42 ) 

A 2 

= - —(a*n*,n*) 

A + 

A = A * 

+ (a*n*,n*) 


that 


(5.43) 
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But it follows from the definition of S v that (y —x*,n*} = —rjcosO, therefore 
we conclude that 


A = — 


\+r)cos9 


(5.44) 


(a*n*,n*)' 

As a conclusion, we have shown that the quadratic approximation of 
the potential assumes its minimum on S v at the point 


X+ycosO , 

y = x* + H+ n*, 

(a*n*,n*) 

(5.45) 

x* satisfies 


iV*(y - x*) = -A \{y - x*). 

(5.46) 


Since TV* is the linearisation of the vector field driving the instanton dynamics, 
we are led to the following remarks. 

(i) The tangent hyperplanes to the equipotential hypersurfaces along the 
instanton trajectory are all parallel to each other (in the neighbour¬ 
hood of the saddle point), and to the separating hyperplane between 
the basins of attraction of X\ and X 2 , see Figure 01 

(ii) The eigenvalue of iV* associated with the direction of the instanton is the 
opposite of the eigenvalue of M* associated with the unstable direction. 


5.4.2. Gaussian approximation of /(. Since y — x* is of order 1 / y/e. we 
may use the approximation 

J ^ ex p( _ —— dS{y), (5.47) 

and (15.461) yields 

(iV*(j/-x*),n*) = -X* + (y- x*,n*) = \\ycos6. (5.48) 

Let us now compute the remaining Gaussian integral. The previous 
paragraph provides the canonical decomposition 

(V - a?*, H+{y - x*)) = (y - y, H*(y - y)) + (y - x*, H+(y - x*)) (5.49) 

for all y € S v . On the one hand, 

/- r, /- ,, A*/? 2 cos 2 9 

(y — x*, HAy — x*)) =—- 7-. (5.50) 

K " (a*n*, n*) V ' 

On the other hand, 


'y^Sr, 


exp 


(y-y,H*(y-y)) 

2e 


dS(y) = / exp 
JveV- 


(v,H+v) 

2e 


d S(y), 
(5.51) 

so that denoting by h £ 1 ) x ( d B the matrix with coefficients (ei,H+ej) 

where (ei,..., e^-i) is an orthonormal basis of V_, we obtain 



(y-y.H^y-y)) ^ 

2e J 


d S(y) 


(27re) 


d -1 


det h 


(5.52) 
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We conclude that 


A +77 cos 9 


(2 


ire) d 1 


det h 


exp 


/ A* ?y 2 cos 2 9 \ 
\ 2 e(a*n*,n*)/ ' 


(5.53) 


5.4.3. Expression of det h. In this paragraph, we prove that det h is connected 
with and det if* through the identity 


det if* = 


Aij_ det h 
{a+n*, n *)' 


(5.54) 


In this purpose, we first establish a variant of the formula (15.101) . Us¬ 
ing (15.81) and the definition of Ai*, we have 


if*Ai* + Ai* T if* + 2if*a*if* = 0. (5.55) 

Multiplying both sides by if* - , we deduce that Af* solves the stationary 
Lyapunov equation 

Ai*ff*- 1 +if*' 1 Ai* + 2a* = 0. (5.56) 


We now complete the orthonormal basis (ei,..., e^-i) of V- introduced 
above by adding the normal vector := n*. In this basis, if*, if * -1 and Ai* 
assume the block decomposition 


if* = 


h 

E r 

E 

7 


if*- X = 


* 

F t 

F 

0 


AL* = 


m 

* 

0 

a ; 


(5.57) 

The cf— 1 first coefficients of the last line of Ai* vanish because of the stability 
of V- by Ai*, while the fact that the last coefficient is At immediately follows 
from the identity (15.401) . Evaluating (15.561) with these block decompositions 
yields 


But by Cramer’s rule, 


A+/3 + (a*n*, 


0 


1 

det if 


n*) = 0 . 

(5.58) 

det h, 

(5.59) 


whence (15,54|) . 


5.5. Conclusion 

Introducing the function 

4(r):= vfcl-U xp K) ds ’ (5 - 60) 

so that, for all y € S v , the commitor function obtained Subsection 15.31 writes 



I cos 2 9 
e(a*n*,n*) 


q e (y) = $ 


(5.61) 
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we compile the results of Subsection 15.41 and obtain the following expression 

C st (x*) 


A e T 


f d /2 


■ exp - 


V(xi,x it ) 


® ~V\ 


I A+ cos 2 9 
e(a*n*,n*) 


(27re) d 1 A^ 


(Kv 2 


cos 2 9 


(5.62) 


x A* ?y cos / ———— -- exp —-- 

y | detU*|(a*7i*,n*) \2e(a*n*,n*) 

Using the asymptotic equivalence 

we observe that the terms depending on rj compensate and therefore lead to 

V(x i,x* 


(5.63) 


AJ 


= U st (x*) exp - 


d—2 


a* 


( 2?r ) 

I det iJ* 


(5.64) 


Using the expression of C st obtained in Section [3] and noting that in the 
non-Gibbsianness correction, the instanton must be integrated for i £ K, we 
conclude that the average transition time, given by the inverse of A| 1 _ >$2 , 
has the expression 


F r e I _ / | detff^l 

Xl ^ X2 eio y det Hess x V(xi,xi) 


exp 


r + OO 


F(p t )dt 


x exp 


f V(x i,x) 


(5.65) 


This formula is very similar to the classical Eyring-Kramers formula for re¬ 
versible diffusions, with the quasipotential playing the role of the potential. 
The eigenvalue A+ can be understood either as the unstable eigenvalue of the 
relaxation dynamics, or as the stable eigenvalue of the instanton dynamics 
around the saddle-point. Finally, the extra term involving the instanton tra¬ 
jectory is related with the prefactor of the quasistationary distribution and 
does not depend on the behaviour of the system around the saddle-point. 
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